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ABSTRACT 

We present a new model for the behavior of scattered time-dependent, asym- 
metric near-UV emission from the nearby ejecta of n Car. Using a 3-D hydrody- 
namical simulation of y Car's binary c o lliding winds, we show that the 3-D binary 



orientation derived by iMadura et al.l ( 120121 ) is capable of explaining the asym- 



metric near-UV variability observed in the Hubble Space Telescope Advanced 
C amera for Su r veys/H igh Resolution Camera (HST ACS/HRC) F220W images 



of ISmith et al.l ( 12004bl ). Models assuming a binary orientation with i m 130° to 
145°, u « 230° to 315°, PA Z « 302° to 327° are consistent with the observed 
F220W near-UV images. We find that the hot binary companion does not signif- 
icantly contribute to the near-UV excess observed in the F220W images. Rather, 
we suggest that a bore-hole effect and the reduction of Fe II optical depths inside 
the wind-wind collision cavity carved in the extended photosphere of the pri- 
mary star lead to the time- dependent directional illumination of circum-binary 
material as the companion moves about in its highly elliptical orbit. 

Subject headings: stars: atmospheres — stars: individual (Eta Carinae) — stars: 
mass-loss — stars: variables: general — supergiants 



Introduction 



Multi- wavelength observations obtained over the past twenty years strongly indicate that 
rj Carin ae is a highly eccentric (e ~ 0.9) colliding wind binary with a 2022.7±1.3 day orb ital 
period (jDaminelilll996l : iPittard fc Corcoran! l2002t ICorcoranll2005t iDamineli et al.ll2008al Jbl). 
With a total luminosity >5x 10 6 L & , dominated by the primary star r)A, a luminous blue 



^ased on HST ACS/HRC observations. 
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2001, 2006 



variable (jDavidson fc Humphreyslll997l ). rj Car's total binary mass is > 120 il^f lHillier et al. 



hereafter HOI, H06). 



One topic that remains the subject of debate is the orientation of r\ Car's orbit. Most 



riastron (Damineli et al. 


19971 


Pittard & Corcoran 2002 


Corcoran 2005: Ioine et al 


2005: 


Hamaeuchi et al. 12007: N 


ielsen et al.ll2007 


Damineli et al 


. 2008b: Okazaki et al.ll2008: 


Gull et al. 


2009. 201 ll: Parkin et al. 


2009; 


Groli et al. 


2010bf). However, others olace /7r on the ne 


ar side 
ad ref- 


of ??a at periastron (e.g. 


Falceta-Goncalves & Abraham z 


>009: Kashi k Soker 2009. a 



erences therein). Settling this debate is crucial as a precise set of orbital parameters is key 
for determining the individual stellar masses. 



Smith et al.l (l2004bl . hereafter S04) attempted to constrain the geometry of r\ Car's orbit 
using asymmetric variability seen in near-ultraviolet (NUV) images of the Homunculus neb- 
ula obtained with the Hubble Space Telescope Advanced Camera for Surveys/High Resolution 
Camera (HST ACS/HRC). Alternating patterns of bright spots and 'shadows' observe d on 
opposite sides of i] Car before and after its 2003.5 spectroscopic event are interpreted by lS04 
as being due to a time- variable, asymmetric NUV radiation field that arises from (1) intrinsic 
NUV emission from t^a's outer wind and (2) UV radiatio n from t?b that preferentially escapes 
in directions away from t?a- In the scenario proposed by IS04l a dark shadow appears on the 
opposite side of t?a near periastron because its dense wind blocks t^b's far-UV radiation over 
a large, time-varying solid angle. Using this interpretation, IS04J suggest that r] Car's orbital 
major axis is nearly perpendicular to the observer's line-of-sight, with t?b on the far side 
of i]a before periastron, on the near side after, and orbiting clockwise on the sky (see their 
figure 2). 



Recently, iMadura et al.l (120121 hereafter M12) tightly constrained, for the first time, 
the 3-D orientation of r\ Car's orbit using a 3-D dynamical model for the broad, spatially- 
extende d [Fe HI] emissio n obser ved by the EST Space Telescope Imaging Spectrograph 



(STIS) (jGull et al.l 120091 ). |M12| find that the observer's line-of-sight has an argument of 
periapsis u ~ 240° to 285°, with the binary orbital axis closely aligned in 3-D with the 
Homunculus polar axis at an inclination i « 130° to 145° and position angle on the sky 
PA^ « 302° to 327°, implying that t/b indeed orbits clockwise on the sky. 

Using a 3-D hydrodynamical model of rj Car's binary colliding winds, we show in this 
letter that the orbital orientation derived by |M12| is consistent with th e asy mmetric NUV 
variability observed in the HST ACS/HRC images of rj Car presented in IS04l The model in 
this letter builds on the earlier work of IS04J . but differs in a key respec t, accounting for the 
wind - wind collision (WWC) cavity created by r?B i n t/a's dense wind (IPittard fe Corcoran 



20021 : lOkazaki et al.l 120081 : IParkin et al.l 120091 : lM12h . This cavity reduces the H and Fe II 
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optical depths in line-of-sight to 77A (Groh 2011 and Groh et al. 2012, in preparation, hereafter 
Gil and G12, respectively ), and causes a bore-hole effect ( IMadura k, Owockill2010t iMadura 



20101 ; IMadura et al.ll201ll ). wherein the WWC cavity allows increased escape of continuum 



radiation from the hotter/deeper layers of rj^s extended wind photosphere at phases near 
periastron. The results of this letter provide insights into how/where NUV light escapes 
the 77 Car binary system and the time-dependent illumination of 77 Car's ejecta in various 
directions. 



2. Observations and Modeling 



We use the same difference images of 77 Car as shown in figure 1 of |S04| . to which we 
refer the reader for further details. We exa mine the F220W fil ter NUV images (probing 
the wavelength region from ~ 1800 - 2600A, ISirianni et al.l 120051 ) where the observed NUV 
excess is greatest. Each image frame shows the result of subtracting the average of all six 
observations from the original image at the indicated date. Here, we focus on the centr al 
±2" region where the observed asymmetric brightness changes are most pronounced (S04J ). 



We zero the phase 6 in 77 Car's 5.5-yea r cyc le to JD = 2,452,819.8 with period 2022.7 days 



( iDamineli et al.ll2008al ). As discussed in IS04L most of the observed variability occurs 0.5" 
to 1" from the central stellar source, thus, <p must be corrected to take into account light 
travel time. The phase delay for regions 1" from 77 Car's central source, adopting D = 2.3 kpc 
( 1Smithll2006l ). is A(p ps —0.007. The revised phases corresponding to the individual difference 
images are therefore = 0.865, 0.925, 0.985, 0.003, 0.031, and 0.061. It is assumed for 
simplicity that phase zero of the spectroscopic cycle (the observations) coincides with phase 
zero of the orbital cycle (periastron). In a highly-eccentric binary system like r\ Car, the tw o 
values are not expected to be shifted by more than a few weeks (I Groh et al.l 1201 Oat |M12| ). 
and so this assumption does not greatly affect our overall conclusions. 

We use a 2-D radiative transfer code and a 3-D hydrodynamical model of 77 Car's collid- 
ing winds to interpret the phase-dependent HST ACS/HRC F220W images. To determine 
the influence of t?b, its low-density wind cavity, and the dense WWC- region wall s on the ob 



served spectrum of 77 Car, we use the 2-D radiative transfer models described in I Groh et al. 



( 12010al ) and G11,G12. A 3-D Smoothed Particle Hydrodynamics (SPH) simulation is used 
to understand the effects of orbital motion on the WWC region formed between 77A and 773, 
and the spatial orientation of the WWC surface on the sk y as a function of phase. The 
3-D SPH simulation in this letter is identical to that used in |M12| . with the exception of the 
size of the computational domain, which here is a factor of ten smaller in order to focus on 
the dynamics of the inner WWC zone. The adopted stellar, wind, orbital, and orientation 
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parameters are in Table [H 



3. Origin of the near-UV Flux in the F220W Images 



Before examining the directional illumination of rj Car's circumstellar ejecta, it is im- 
portant to understand the physical origin of the observed NUV flux. According to the 2-D 
radiative transfer model of G11,G12, regions in line-of-sight to the wind cavity carved by t)b 
are exposed to a larger NUV flux (red line of Figure [1]) than regions in line-of-sight to only 
the extended wind of 7?a (blue line of Figure [T]). Between 1800 and 2600A, these modified 
wind models of 7?a have a flux that is roughly an order of magnitude higher than the flux 
of t]b- We assume for r/B a temperature of 35,000 K and luminosity o f 10 6 L Q (H06). Since 
10 6 L & is an upper limit for the luminosity of 77b (IMehner et al.ll2010l ). the flux contribution 
from t]b could very likely be even less than what is shown in Figure [TJ Therefore, we find that 
t/b does not contribute a substantial fraction of the observed F220W NUV flux in directions 
away from rjx and cannot explain the observed F220W NUV excess. Rather, we suggest 
that the NUV excess in the F220W images is regulated by the time-dependent nature of 
the WWC cavity carved by t/b, which has two very important effects on tja's extended wind 
photosphere. 

First, lines-of-sight through the low-density WWC cavity have significantly reduced Fe II 
optical depths (G11,G12). This redu ction is especially pronounced in the F220W spectral 
region, which is full of Fe II lines (IHOq ; G12). Second, t/a's photosph ere is extremely extended 



at U V wavelengths f|H06l) . leading to a significant bore- hole effect ( Madura fe Owockil 12010 



Madura! l2010t iMadura et al.l 1201 ll ) at phases around periastron, wherein the WWC cavity 
creates a hole, allowing increased escape of continuum radiation from the exposed hotter, 
deeper layers of t/a's photosphere. 

Because r/ B is located well within the WWC cavity, <10% of the F220W flux that 
reaches circumstellar ejecta exposed to the cavity comes from r/B- Therefore, the time- 
dependent directional exposure of circumstellar material to excess F220W NUV flux most 
likely depends on the phase- dependent spatial orientation of the WWC cavity, and not just 
the obstruction of t/b's NUV flux by r}\s wind. 



4. The Lighthouse Effect: Constraints on the Orbital Orientation Parameters 



Using our 3-D SPH model we investigated the phase-dependent spatial orientation of 
the WWC cavity for different orbital orientations with the goal of determining which ori- 
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Fig. 1. — Comparison of model spectra in the range 1800 — 2600A, corresponding to the 
wavelength region probed by the HST ACS/HRC F220W NUV images. All spectra have 
been smoothed to a resolution of R = 600 to make differences between models m ore c learly 
visible. Green Line: Spectrum of t/b computed using the 1 -D CMFGEN model of lH06l . Blue 
Line: Unmodified model spectrum of rjx from iHOll ; lH06l . Red Line: 2-D model spectrum 
from G12 at apastron, which includes the WWC cavity and its density-enhanced walls. 
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entations are capable of explaining the HST F220W NUV images of |S04| . For simplicity, 
each orientation discussed in this paper assumes that the orbital axis of the binary system 
is closely aligned in 3-D with the polar axis o f the Homunculus nebula at an inclination 
i = 138° and position angle on the sky of 312° ( 1M12I ). 



A quantitative description of changes in the amplitude of the flux in the F220W images 
is deferred to future work as a 3-D radiative transfer code is necessary for such an analysis. 
Below, we compare the spatial location of features in the individual F220W images to those 
expected from the model based on the 3-D orientation of the WWC cavity opening. Compass 
directions (NE = north-east, etc.) are used in the descriptions below. 



Figures [2] and [3] show that the binary orientation proposed by lM12l is capable of reason- 
ably explaining the observed time-variable NUV 'excess' emission seen in the F220W images. 
Hereafter, the terms 'NUV excess' and 'NUV deficit' refer to an increase and decrease, re- 
spectively, in the F220W NUV flux as compared to the average NUV flux of phases between 
0.865 and 0.061. 

At (j) = 0.865 and 0.925, the observed NUV excess emission in the F220W images is 
brig htest in regions in line-of-sight to the central stellar source and is elongated from NE to 
SW (JSOJ). t]b is far enough from t)a at these phases that orbital velocities are low and the 
WWC cavity maintains a simple axisymmetric cone - like shape (top two rows of Figure 2; also 
Okazaki et alll2008[ IParkin et allhood l201l[ lM12h . For u « 260°, the WWC cavity is open 
mostly toward the observer and is elongated in directions from the NE to the SW on the sky, 
causing the Fe II optical depths in those directions to be reduced (G12). As a result, material 
in these directions should be exposed to a higher F220W NUV flux than the one- year average 
flux of phases between 0.865 and 0.061. The NUV excess thus occurs in these directions at 
these two phases because at later times the inner WWC cavity has a different 3-D spatial 
orientation (Figure ED- Moreover, the models predict that because orbital velocities are low, 
the orientation of the WWC cavity should not change much between <fi = 0.865 and 0.925, 
implying that the spatial orientation and amount of observed NUV excess should also not 
change much. The F220W images show that this is the case. 

By (f) = 0.985 (bottom row of Figure[2]), the spatial orientation of the inner WWC cavity 
has changed, pointing more in directions to the SW on the sky. The WWC cavity at this 
phase is warped due to orbital motion. The F220W image shows an observed NUV excess 
i n di rections to the NW, SW, and SE, with darker, below-average-flux regions to the NE 
(1S04J ) . According to our model and interpretation, material to the NW, SW, and SE should 
have a NUV excess as it is exposed to a heavily modified primary wind. The darker region 
to the NE is likely due to there being more primary wind material in this direction at this 
phase compared to the average of all phases. 
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Fig. 2. — Comparison of observed HST F220W images (left column, IS041 ) with the 3-D SPH 
model (right two columns) discussed in the text at phases (rows, top to bottom) (f> = 0.865, 
0.925, and 0.985. Middle column: 3-D isovolume renderings of the SPH simulation assuming 
the parameters in Table [TJ illustrating the orientation of the WWC cavity (dark purple) 
carved in rj^s extended wind (in red) as it would appear projected on the sky. The surface 
of the WWC cavity is color coded to radius, i.e. dark purple indicates material at larger 
radii (~ 75 to 125 AU) from the central stars, while green indicates material near the apex 
of the WWC region (< 15 AU from r]A for these phases). Right column: Same as middle 
column, but with material below the orbital plane removed in order to show the complex 
dynamics of the inner WWC cavity. The inset in the lower left shows the orientation of the 
binary orbit (yellow) projected on the sky, along with the semi-major (x, red), semi-minor 
(y, green), and orbital (z, blue) axes. The black arrows indicate the clockwise orbital motion 
of the stars. North is up and east is left in all panels. 
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Fig. 3. — Same as Figure El but for phases (top to bottom) <fi = 0.003, 0.031, and 0.061. 
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At = 0.003, the inner WWC region has started to take on a spiral shape and points 
mainly in directions to the S on the sky and partly away from the observer (top row of 
Figure [3D- Interestingly, the F220W image shows NU V excesses to the S and E on the sky, 



with dark NUV-deficit regions to the N and W (lS04j ) . At = 0.003, there is a significant 
bore-hole effect concentrated in directions to the S on the sky. The primary wind is carved 
to the S as well, which should allow the increased flux from the bore-hole to reach material 
in this direction. Since there is no bore-hole effect to the S at the other phases, we expect a 
NUV excess in material there at = 0.003. Similar reasoning explains why there is also a 
NUV deficit to the N and W. 

At = 0.0 31 (middle row of Figure [3]), the inner ±2" region of the F220W image is 



very dark ((SOJ). The inner WWC cavity gets highly distorted and the additional NUV 



radiation that wo uld es cape from the inner layers of //a becomes embedded and trapped 



by its dense wind (IM12I) . Therefore, a NUV deficit and a large dark region appears in the 



F220W difference image. 

By = 0.061 (bottom row of Figure [3]), orbital speeds have decreased and the WWC 
cavity has started to increase in size in directions to the NNE on the sky, significantly carving 
the wind of t?a in this direction. Consequently, one would expect a significant decrease in 
the amount of Fe II absorption in directions to the NNE compared to the average, and thus 
an observed NUV excess. There is also more wind material from r]\ in directions to the S 
compared to the average, which should cause a NUV deficit. The F220W image show s that 



there is indeed an observed NUV excess to the NNE and a dark region to the S ((SOJ) 



We find that only binary orientations with 230° < u < 15° are capable of reasonably 
explaining the extended (~ 0.2" to 2") NUV excess emission observed in the F220W images. 
Other orientations appear to have extreme difficulty explaining the NUV variability. For 
example, at an orientation of u = 90° (top row of Figure 0J, during most of the binary 
orbit, the WWC cavity is open and pointing away from the observer and to the E/SE. Thus, 
between = 0.865 and 0.985, the observer should detect the NUV radiation field from the 
unmodified wind of i]a on the central stellar source, and a possible NUV excess to the E/SE. 
At = 0.061, due to the wrapping of the WWC region, a NUV excess would be expected 
to the S. Yet, the spatial location of the observed bright and dark spots is nearly the exact 
opposite. 

However, we find that the NUV variability on the central stellar source (central ±0.2") 
additionally constrains the orbital orientation to values of u ~ 230° to 315°. The F220W 
images at = 0.865, 0.925, and 0.985 show an observed NUV excess o n the central stellar 
source, while = 0.003, 0.031, and 0.061 show a deficit (|S04 ; see also Martin et al. 2006 



Mehner et al .1120 111 ). In order to have a NUV excess on the central source before periastron, 
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4p = 0.061 




RA Offset (arcsec) 



Fig. 4. — Similar to Figures [2] and El but for orbital orientations assuming to = 90° (top row) 
and co = 0° (bottom two rows) at selected phases. The red wind of t]a in the renderings of 
the last two columns has been made semi-transparent in order to allow the reader to see the 
dark purple surface of the WWC region, which opens away from the observer for these two 
orientations. 



-li- 



the WWC cavity should be open mostly toward the observer at these phases. NUV deficits 
after periastron are due to increased amounts of primary wind material (compared to the 
average) at these times. This appears to be the case for u ~ 230° to 315°. In contrast, 
assuming u = 0°, for example, the WWC cavity opens away from the observer at phases 
= 0.865 to 0.003 (Figure H]). At (j) = 0.061, the WWC spirals in between the observer and 
?7a (bottom row Figure @J. This would lead one to expect a NUV deficit before periastron 
and an excess after if uj = 0°, which is not seen in the observations. 

Moreover, binary orientations that lie significantly outside the range u ~ 230° to 315° 
have great difficulty explaining other multi-wavelength diagnostics of rj Car since such ori- 
entations place th e observer's line-of-sight through tm's optically-thick wind for most of 



Groh et al. 


2010b 


proposed by 


M12| 



M12 



Pittard &: Corcoran 2002; Okazaki et al 



2008: Parkin et al. 2009, 2011 



G11,G12). In contrast, the 3-D orientation and direction of orbit 



those analyzed here. 



The results of this letter go well beyond constraining the orientation of rj Car's binary 
orbit. It is clear that the motion of the WWC cavity as t?b moves about in its highly 
elliptical orbit leads to an important 'lighthouse effect' in rj Car, wherein circum-binary 
ejecta is exposed to a time- and direction-dependent modified NUV radiation field of t]a 
caused by a bore-hole effect and decrease in Fe II optical depths. This lighthouse effect is 
crucial for understanding how NUV light escapes the rj Car system and the phase-dependent 
illumination of distant ejecta in different directions. 

The lighthouse effect also provides a valuable diagnostic f or he lping constrain the exact 
timing of periastron. Based on the available observations of IS04j and the simple model in 
this letter, periastron should occur between (ft = 0.985 and (p = 0.031, most likely very close 
to (j) — 0.003 since this is when some NUV excess is still visible to the SE in the difference 
images (i.e. before distortion of the WWC cavity causes the NUV radiation to be trapped 
by t^a's dense wind). 

Future spatially-resolved observations with better time sampling, together with detailed 
3-D radiative transfer models, can help place much tighter constraints on the exact timing of 
periastron, and possibly the orbital eccentricity. We emphasize that future observations of 
77 Car should focus not only on phases around = 0, but also on the extended recovery period 
up until 4> ~ 0.2, during which time the WWC cavity is increasing in size and reestablishing 
its axisymmetric cone-like shape. Future monitoring of these phases is crucial for determining 
when various forms of radiation can escape via the WWC cavity in directions away from rj\. 
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Table 1: Adopted Model Stellar, Wind, Orbital, and Orientation Parameters 



Parameter t/a ??b 

Mass (M ) 90 30~ 

Mass-Loss Rate (M yr" 1 ) 10~ 3 KT 5 

Wind Terminal Velocity (km s^ 1 ) 500 3000 

Orbital Period (days) 2024 

Orbital Eccentricity e 0.9 

Semi-major Axis Length a (AU) 15.4 

Orbital Inclination i 138° 

Argument of Periapsis u 260° 

Position Angle on Sky of Orbital Axis, PA Z 312° 



